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ABSTRACT 

We consider astrophysically relevant nonlinear MHD dynamo at large Reynolds numbers (Re). 
We argue that it is universal in a sense that magnetic energy grows at a rate which is a constant 
fraction C'e of the total turbulent dissipation rate. On the basis of locality bounds we claim that 
this "efficiency of small-scale dynamo", Ce, is a true constant for large Re and is determined only 
by strongly nonlinear dynamics at the equipartition scale. We measured Ce in numerical simulations 
and observed a value around 0.05 in highest resolution simulations. We address the issue of Ce being 
small, unlike Kolmogorov constant which is of order unity. 
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1. INTRODUCTION 

MHD turbulence is ubiquitous in astrophysical and 
space environments (Goldstein et al 1995; Armstrong et 
al 1995). Reynolds numbers are, typically, very high, 
owing to large scales of astrophysical processes com- 
pared to small dissipative scales. One of the central 
questions of MHD dynamics is how initially unmagne- 
tized well-conductive fluid generates its own magnetic 
field, known broadly as "dynamo" . Turbulent dynamo 
has been roughly subdivided into "large-scale dynamo" 
and "small-scale" or "fluctuation dynamo" depending on 
whether magnetic fields are amplified on scales larger 
or smaller than turbulence outer scales. Although a 
few "no-dynamo" theorems has been proved for flows 
with symmetries, a generic turbulent flow, which pos- 
sess no exact symmetries was expected to amplify mag- 
netic fields by stretching, due to particle separation in 
a turbulent flow. For the large-scale field, a "twist- 
stretch-fold" mechanism was introduced (Vainshtein & 
Zeldovich 1972). Turbulent flow possessing perfect sta- 
tistical isotropy can not generate large-scale magnetic 
field. Observed large-scale fields, such as in disk galax- 
ies, are generated when statistical symmetries of turbu- 
lence are broken by large-scale asymmetries of the sys- 
tem, such as stratification, rotation and shear (see, e.g., 
Vishniac & Cho 2001; Kapyla et al 2009). Since these 
symmetries are only weakly broken, large-scale dynamo 
is slow. Small-scale dynamo does not suffer from this re- 
striction and can be fast. Kinematic small-scale dynamo, 
which ignores the backreaction of the magnetic field has 
been studied extensively before (Kazantsev 1968; Kul- 
srud & Anderson 1992). However, kinematic dynamo 
is irrelevant for astrophysical environments, because it 
grows exponentially and becomes inapplicable at small- 
est turbulent timescales which are tiny by astrophysical 
standards. Additionally, kinematic dynamo have posi- 
tive spectral index (typically 3/2) at all scales, which is 
incompatible with observations in galaxy clusters (Laing 
et al 2008) which clearly indicate steep spectrum with 
negative power index at small scales. Due to preexisting 
astrophysical fields small-scale dynamo starts in nonlin- 
ear regime. It was discovered numerically that small- 
scale dynamo continues to grow after kinematic stage, 
producing steep spectrum at small scales and signifi- 



cant outer-scale fields (Haugen et al 2004; Brandenburg 
& Subramanian 2005; Cho et al 2009; Ryu et al 2008). 
Furthermore, MHD turbulence produces turbulent diffu- 
sivity (aka "/3-effect"), which is essential for large-scale 
dynamo (Kapyla et al 2009) and reconnection (Lazarian 
& Vishniac 1999; Eyink 2010). Saturation of small-scale 
dynamo seems to be independent on Re and Pr as long as 
Re is large (Haugen et al 2004) and the magnetic energy 
growth rate could be constant (Schekochihin & Cowley 
2007; Cho et al 2009; Beresnyak et al 2009; Ryu et al 
2008). Small-scale dynamo is faster than large-scale dy- 
namo in most astrophysical environments and magnetic 
energy grows quickly to equipartition with kinetic mo- 
tions, with the largest scales of such field being a fraction 
of the outer scale of turbulence. Subsequently, these tur- 
bulent fields are slowly ordered by mean-field dynamo, 
with turbulent diffusivity of MHD turbulence playing es- 
sential role. In this paper we provide sufficient analytical 
and numerical argumentation behind the universality of 
nonlinear small-scale dynamo. 

2. NONLINEAR SMALL-SCALE DYNAMO 

We assume that the spectra of magnetic and kinetic 
energy at a particular moment of time are similar to 
what is presented on Fig. 1. Magnetic and kinetic spec- 
tra cross at some "equipartition" scale 1/fc*, below which 
both spectra are steep due to MHD cascade (see, e.g., 
Goldrcich & Sridhar 1995; Beresnyak 2011). This is sug- 
gested by both numerical evidence (Beresnyak & Lazar- 
ian 2009b; Cho et al 2009) observations of magnetic fields 
in clusters of galaxies (Laing et al 2008). At larger scales 
magnetic spectrum is shallow, k a , a > 0, while kinetic 
spectrum is steep due to a hydro cascade. Most of the 
magnetic energy is concentrated at scale 1/k*. We desig- 
nate Ck and Cm as Kolmogorov constants of hydro and 
MHD respectively. The hydrodynamic cascade rate is e 
and the MHD cascade rate as t^- Due to the conserva- 
tion of energy in the inertial range, magnetic energy will 
grow at a rate e — e^- We will designate Ce — (e— £2) A as 
an "efficiency of small-scale dynamo" and we will argue 
that this is a true constant, since: a) turbulent dynam- 
ics is local in scale in the inertial range; b) neither ideal 
MHD nor Euler equations contain any scale explicitly. 
Magnetic energy, therefore, grows linearly with time if 
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Fig. 1. — A cartoon of kinetic and magnetic spectra in small- 
scale dynamo, at a particular moment of time when equipartition 
wavenumbcr is k* . 

€ = const. The equipartition scale 1/fc* will grow with 
time as i 3 / 2 (Beresnyak et al 2009). This is equivalent 
to saying that small-scale dynamo saturates at several 
dynamical times at scale 1/fc* and proceeds to a twice 
larger scale (Schekochihin & Cowley 2007). If magnetic 
energy grows approximately till equipartition (Haugen 
et al 2004; Cho et al 2009), the whole process will take 
around several dynamical timescales of the system, or 

more quantitatively, (C 3 ^ 2 /Ce)(L/v l ). 

3. LOCALITY OF SMALL-SCALE DYNAMO 

We will use "smooth filtering" approach with dyadic- 
wide filter in k-space (Aluie & Eyink 2010; Eyink 2005). 
We designate a filtered vector quantity as where k 
is a center of a dyadic Fourier filter in the range of wave 
numbers [fc/2,2fc]. The actual logarithmic width of this 
filter is irrelevant to further argumentation, as long as it 
is not very small. We will assume that the vector field a 
is Holger-continuous with some exponent and designate 
cifc = (|al fc ] I 3 ) 1 / 3 which has to scale as fc CT3 , e.g., fc _1//3 
for velocity in Kolmogorov turbulence. The energy cas- 
cade rate is e = C K 3 ^ 2 kvl, where here we defined Kol- 
mogorov constant Ck by third order, rather than second 
order quantities. We will keep this designation, assuming 
that traditional Kolmogorov constant could be used in- 
stead. We use spectral shell energy transfer functions 
such as T vv (p,k) = -( V [ fe ](v • V)vW ),T W+W+ (p,k) = 

-(w+ lfe] (w- • V)w+ [pl ) (Alexakis et al 2005), applica- 
ble to incompressible ideal MHD equations, where 
are Elsasser variables and v, b and w" 11 are measured in 
the same Alfvcnic units. We will use central frequency 
k and study "infrared" (IR) transfers from p <C fc, and 
"ultraviolet" (UV) transfers, from 9 > fc. We will pro- 
vide bounds on |T|, in units of energy rate as in Aluie & 
Eyink (2010); Eyink (2005) and relative volume-averaged 
bounds which are divided the actual energy rate and are 
dimensionless. We will consider three main wavevector 
intervals presented on Fig. 1: fc <C fc* ("hydro cascade"), 
fc - fc* ("dynamo") and fc > fc* ("MHD cascade"). 

4. MHD CASCADE, K 2> K*. 

The only energy cascades here are Elsasser cascades 
and, by design of our problem, w + and w~ have the same 
statistics, so we will drop ±. For an exchange withp <C fc 
band, for T^^l, using Holger inequality and wavenum- 
ber conservation we get an upper bound of pw v w\ and 
for q ^> fc band it is kw q Wk, these bounds are asymptoti- 



TABLE 1 
Transfers and upper limits 



Transfers 




p < k 


q > k 


T vv (p,k) = 


-<vW(v . v)vW) 


pvpv'i 


kv k v' 2 


T bb (p,k) = 


-(bW (v ■ V)bW> 


PVpVkbk 


kbkVqbq 


T vb (p,k) = 


(bW(b.V)vW) 


pv p b\ 


kbkVqbq 


T bv (p,k) = 


(vW(b-V)bW) 


pb P v k b k 


kv k b 2 q 


T w+w +(p,k) = 


-(w+ W (w- • V)w+ 


p] ) pw v w\ 


kw k w 2 



cally small, see Eyink (2005). For the full list of transfers 
and limits refer to Table 1. The relative bound should 
be taken with respect to C M kwf, where Cm is a Kol- 
mogorov constant for MHD, from which we get that most 
of the energy transfer with the [fc] band should come from 

[fcC~ 9/4 ,fcC^ 4 ] band, see Beresnyak (2011). The global 
transfers between kinetic and magnetic energy must av- 
erage out in this regime, nevertheless, the pointwise IR 
and UV transfers can be bounded by pb p Vkbk and kb 2 q Vk 
and are small (Eyink 2005). 

5. HYDRO CASCADE, K <C K* 

Despite having some magnetic energy at these scales, 
most of the energy transfer is dominated by velocity field. 
Indeed, \T VV \ is bounded by pv p vf for p -C fc and by 
kv^Vk for q fc. Compared to these, \Tb v \ transfers are 
negligible: pbpVkbk and kb 2 q Vk- For magnetic energy in 
p <C fc case we have \T v b\ and |Xf,fc| transfers bounded 
by pvpb 2 ,, pbpVkbk and for g > fc case \T v b\ and \Tbb\ 
are bounded by kbkv q b q . Out of these three expressions 
the first two go to zero, while the third goes to zero if 
a — 2/3 < or have a maximum at q — fc* if a — 2/3 > 0. 
This means that for the transfer to magnetic energy we 
have IR locality, but not necessarily UV locality. Note 
that magnetic energy for fc <C fc* is small compared to 
the total, which is dominated by fc — fc*. We will assume 
that a — 2/3 > and that the spectrum of bk for fc < fc* 
is formed by nonlocal \T vb \ and \T b b\ transfers from fc*, 
namely magnetic structures at fc are formed by stretching 
of magnetic field at fc* by velocity field at fc. Magnetic 
spectrum before fc* is, therefore, nonlocal and might not 
be a power-law. 

6. DYNAMO CASCADE K = K* 

In this transitional regime our estimates of Elsasser 
UV transfer and kinetic IR transfer from two previous 
sections will hold. We are interested how these two are 
coupled together and produce observable magnetic en- 
ergy growth or decay. IR p <C fc* \T v b\ and |Tf,b| transfers 
will be bounded by pv v b\, and pb v Vk* bk* , which go to 
zero, so there is a good IR locality. Ultraviolet transfers 
will be bounded by k*bk*b q v q . This quantity also goes 
to zero as q increases, so there is an UV locality for this 
regime as well. Let us come up with bounds of relative lo- 
cality. Indeed, the actual growth of magnetic energy was 
defined as = e — e 2 = C E C K 3/2 kvl So, p « fc* IR 

bound is k*C 3 E /2 C K 9/4 and UV bound is k*C ? 3/2 C% A . 
We conclude that most of the interaction which result 
in magnetic energy growth must reside in the wavevec- 
tor interval of k*[C 3 J 2 C K 9/i , C E Zj2 C%\ Numerically, 
if we substitute Ck — 1-6, Cm — 4.2, Ce — 0.05 we get 
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TABLE 2 

Three-dimensional MHD simulations 



Run 


n 


N 3 


Dissipation 


(e) 


Re 


C E 


Ml-6 


6 


256 3 


-7.6 ■ W~ 4 k 2 


0.091 


1000 


0.031 ± 0.002 


M7-9 


3 


512 3 


-3.0 ■ 10- 4 k 2 


0.091 


2600 


0.034 ± 0.004 


M10-12 


3 


1024 3 


-1.2 ■ 10- 4 fc 2 


0.091 


6600 


0.041 ± 0.005 


M13 


1 


1024 3 


-1.6 ■ 10- 9 fc 4 


0.182 




0.05 ± 0.005 


M14 


1 


1536 3 


-1.5 • 10~ 15 fc 6 


0.24 




0.05 ± 0.005 



the interval of k* [0.004, 2000]. So, despite being asymp- 
totically local, small-scale dynamo can be fairly nonlocal 
from practical standpoint. 

To summarize, kinetic cascade at large scales and MHD 
cascade at small scales in the inertial range are domi- 
nated by local interactions. The transition between ki- 
netic cascade and MHD cascade is also dominated by 
local interactions, and since ideal MHD equations do not 
contain any scale explicitly, the efficiency of small-scale 
dynamo Ce is a true universal constant. Note that Ce 
relates energy fluxes, not energies, so this claim is unaf- 
fected by the presence of intermittency. Magnetic spec- 
trum at k <C k* is dominated by nonlocal triads that 
reprocess magnetic energy from k — k* and, since this 
part of the spectrum contains negligible magnetic energy, 
our universality claim is unaffected by this nonlocality. 

7. NUMERICAL RESULTS 

We performed numerical simulations of statistically ho- 
mogeneous isotropic small-scale dynamo by solving MHD 
equations with stochastic non-helical driving and explicit 
dissipation with Pr m = 1. The details of the code and 
driving are described in detail in our earlier publications 
(Beresnyak 201 f; Beresnyak & Lazarian 2009) and Ta- 
ble f shows simulation parameters. We started each sim- 
ulation from previously well-evolved driven hydro simu- 
lation by seeding low level white noise magnetic field. 
We ran several statistically independent simulations in 
each group and obtained growth rates and errors from 
statistical (or sample) averages. In all simulations, ex- 
cept MI4, the energy injection rate was controlled. Fig 2 
shows sample-averaged time evolution of magnetic en- 
ergy. Growth is initially exponential and smoothly tran- 
sition into the linear stage. Note, that scatter is initially 
small, but grows with time, which is consistent with the 
picture of magnetic field growing at progressively larger 
scales and having progressively less independent realiza- 
tions in a single datacube. 

8. EFFICIENCY OF SMALL-SCALE DYNAMO 

Our C'e is much smaller than unity. One would ex- 
pect a quantity of order unity because this is a univer- 
sal number, determined only by strong interaction on 
cquipartition scale. If wc refer to the ideal incompressible 
MHD equations, written in terms of Elsasser variables, 
c^w* + S(w T • V)w ± = 0, dynamo could be understood 
as decorrelation of w ± which are originally perfectly cor- 
related with w + = w~ in the hydrodynamic cascade. In 
our case this decorrelation is happening at the cquipar- 
tition scale k* . Being timc-dcpcndcnt, this decorrelation 
propagates upscale, while ordinary energy cascade goes 
downscale. The small value of Ce might be due to this. 
A different viewpoint on small-scale dynamo is a picture 
when magnetic field is being stretched by kinetic motions 
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Fig. 2. — Magnetic energy growth observed in simulations Ml-6, 
M7-9 and M10-13. We used sample averages which greatly reduced 
fluctuations and allowed us to measure Ce with sufficient precision. 

and, therefore, amplified. With small or negligible diffu- 
sion the field becomes extremely tangled (Schekochihin 
et al 2004), quite contrary to the observations and nu- 
merical experiments, where magnetic field has an energy- 
containing scale which is larger than the scale of resistive 
dissipation. In the nonlinear regime, therefore, the in- 
tuitive picture with stretching with resistive diffusion is 
completely wrong. It is inevitable to appeal to turbulent 
diffusion which helps to create larger-scale field. Both 
stretching and diffusion depend on turbulence proper- 
ties at the same designated scale 1/fc*, therefore in the 
asymptotic regime of large Re one of these processes must 
dominate. As we know Ce is small so stretching and dif- 
fusion are close to canceling each other. 

9. KINEMATIC DYNAMO RATES 

A better studied and understood kinematic dynamo 
might shed some light on the problem of small Ce- In 
the kinematic regime, when we neglect Lorentz force in 
the MHD equation, the growth is exponential and the 
rate is expected to come from fastest shearing rates of 
smallest turbulent eddies. Observed rates, however, arc 
smaller which was interpreted as competition between 
stretching and turbulent mixing (Eyink et al 2011). In 
our simulations, in the kinematic regime of M7-9 we ob- 
served growth rate 77^ = 0.0326, where t v = {v/e) 1 / 2 is 
a Kolmogorov timescale, which is consistent with (Hau- 
gen et al 2004; Schekochihin ct al 2004). With minimum 
timescale, T m - ln « 9t Vi 7T m i n = 0.3, which is still small. 

Kazantsev-Kraichnan model (Kazantsev 1968; Kulsrud 
& Anderson 1992) predicts 7T m i n ~ 1. This model, how- 
ever, uses ad-hoc delta-correlated velocity which does not 
correspond to any dynamic turbulence and its statistics is 
time-reversible as opposed to time-irreversible real turbu- 
lence. Time irreversibility of hydro turbulence actually 
mandates that fluid particles separate faster backwards 
in time, since (vjjj) = — 4/5el is negative. 

In order to study the interplay of stretching and diffu- 
sion, we performed several simulations of kinematic dy- 
namo forward and backward in time. We followed full 
three-dimensional evolution of v and b and approximated 
"backward in time" by reversing velocity direction. Ini- 
tial condition for magnetic field was typically random 
noise. Since we couldn't reverse viscous losses in DNS, 
we used viscosity v = 0, but magnetic diffusivity rj > 0. 
In the first set of simulations we set initial velocity as v 
and — v from evolved viscous runs. The growth rates is 
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Fig. 3. — Evolution of magnetic energy forward and backward 
in time (dashed is with reversed x-axis). Inset: a naive simulation 
with initial state +v and — v from Navier-Stokes simulation (not 
a rigorous backward in time simulation). 

shown on inset of Fig. 3. Quite surprisingly, the "back- 
ward" simulation did not produce any growth for sev- 
eral dynamical times. Unexpectedly, simply reversing 
velocity has such a profound impact on kinematic dy- 
namo, despite spectra being very close to each other, 
suggesting that it is not only the spectrum that deter- 
mines growth but rather the actual statistical properties 
of velocity, which will determine whether stretching or 
diffusion wins, i.e. if there is a dynamo or a no-dynamo, 
even in the simple kinematic case. In this simulation 
we observed a typical k 2 "thermal pool" at the end of 
velocity spectrum which had shortest timescales. "Ther- 
mal pool" was clearly time- irreversible, unlike the true 
physical thermal pool, consistent with (Ray et al 2011). 

The next series of simulations were reproducing an ac- 
tual backward in time dynamics. In order to achieve 
this we evolved initial state for a fairly short time with 
v = and then we evolved it for the same time revert- 
ing velocity with v = and confirmed that final state is 
close to initial state, due to reversibility of truncated Eu- 
ler equations. The results for dynamo growth is shown 
of Fig. 3. We see that backward dynamo is faster by a 
factor of 2.0 ± 0.1, which is actually consistent with the 
ratio of particle diffusion forward and backward in time 
in (Eyink 2010). This result again reinforces our state- 
ment that dynamo is a result of competing mechanisms 



of turbulent stretching and turbulent diffusion and the 
outcome depend on statistics of velocity other than just 
velocity spectrum. 

10. DISCUSSION 

Linear growth regime that we considered in this Let- 
ter provides a fastest growth rate for magnetic energy, 
Ce£, which is actually just a constant fraction of abso- 
lute energetic limit of e. This linear growth starts with 
t = in astrophysical enviroments with very high Re, as 
long as turbulence is well-developed. Exponential growth 
regime, that is relevant for finite Re and small initial seed 
fields provide much slower growth. Combined with finite 
time that is available in simulated astrophysical objects 
this may result in severe underestimation of the magnetic 
energy and, therefore, of dynamical importance of mag- 
netic fields. For example, cosmological simulations of 
halo collapse, that have rather modest Re, could be un- 
derestimating magnetic energy (Sur ct al. 2010). A res- 
olution study must confirm convergence of magnetic en- 
ergy with resolution. If exponential growth is observed, 
magnetic energy is underestimated. 

In this Letter we discussed the issue of Ce being small 
and interpreted this as a competition between turbulent 
stretching that enhances magnetic field and turbulent 
diffusion that tries to dissipate it. Indeed, the nega- 
tive values of Ce, from theoretical viewpoint, are also 
allowed. This would be a situation of non-dynamo in the 
nonlinear regime. Isotropic homogeneous Kolmogorov 
turbulence does produce a growth and most astrophysical 
enviroments that are well-magnetized. However, turbu- 
lent magnetic diffusion can dominate in situations when 
turbulence is very different from isotropic and homoge- 
neous and have finite Re. This might explain why dy- 
namo experiments sometimes generate larger fields in the 
absense of turbulence (Colgate et al. 2011). 
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Gregory Falkovich, Alex Lazarian, Hui Li and Hao Xu for 
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lowship and Humboldt Fellowship. Computations were 
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TG-AST080005N. 
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